N = [2:65,70:5:200]';

x = load('tableData.dat');

Hnumerical = x(:,2) ./ 2;

o = optimset('MaxFunEvals',4000,'MaxIter',4000,'Display','off');

x1 = lsqcurvefit(@fitFunctionForH1415,[0,0,0,0,0,0],[2:65]',Hnumerical(1:64),[],[],o);
x2 = lsqcurvefit(@fitFunctionForH,[0,0,0,0,0,0],[2:65]',Hnumerical(1:64),[],[],o);


plot(N,fitFunctionForH1415(x1,N),N,fitFunctionForH(x2,N),N,Hnumerical);
legend('14/15','1','numerical');